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Abstract 

The quantum dynamics of one-dimensional bosonic Josephson junctions with attractive and 
repulsive interparticle interactions is studied using the Bose-Hubbard model and by numerically- 
exact computations of the full many-body Hamiltonian. A symmetry present in the Bose-Hubbard 
Hamiltonian dictates an equivalence between the evolution in time of attractive and repulsive 
Josephson junctions with attractive and repulsive interactions of equal magnitude. The full many- 
body Hamiltonian does not possess this symmetry and consequently the dynamics of the attractive 
and repulsive junctions are different. 
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Quantum dynamics of interacting Bose-Einstein condensates is an active and lively 
research field LU, Here, one of the basic problems studied is the dynamics of 

tunneling of interacting Bose-Einstein condensates in double-wells, which in this con- 
text are referred to as bosonic Josephson junctions. The dynamics of bosonic Joseph- 



son junctions has drawn much attention bot 



Refs. 5, 6, 3, S S 10, 11, 12, Q. 
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1 theoretically and experimentally, see, e.g., 
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18| and references therein. 



In this work we would like to compare the dynamics of one-dimensional bosonic Josephson 
junctions with attractive and repulsive interparticle interactions. Explicitly, we will compare 
systems where the magnitude of attractive and repulsive interactions is alike. We prepare 
the interacting bosons in one well, and then monitor the evolution of the systems in time. 
We compute and analyze the respective dynamics both within the two-site Bose-Hubbard 
model often employed for this problem as well as within the full Hamiltonian of the systems. 
The main result of this work, shown both analytically and numerically, is that within the 
Bose-Hubbard model the dynamics of the attractive and repulsive junctions is equivalent. In 
contrast, the dynamics of attractive and repulsive junctions computed from the full many- 
body Hamiltonian are different from one another. As a complementary result we provide 
here for the first time in literature numerically-exact quantum dynamics of an attractive 
Josephson junction, thus matching our recent calculations on repulsive Josephson junctions 

m 

We begin with the two-site Bose-Hubbard Hamiltonian: 

h(U) = -J (i[b R + P R b L ) + | (b{b{b L b L + b R b R b R b R ) . (1) 

Here and hereafter all quantities are dimensionless. We remind the reader that the 
Bose-Hubbard Hamiltonian ([T]) is derived from the full many-body Hamiltonian 7i = 
J dx& '(x)h(x)4f '(x) + ^ J dx& (x)& (x)4f (x)4f (x) , where h(x) is the one-body Hamil- 
tonian and Ao the interparticle interaction strength. Explicitly, by restricting the field 
operator to a sum of two terms 4?(x) = b L <p L (x) + b R <fi R (x), where and 4> R (x) are 

the left- and right-localized Wannier functions, substituting into the many-body Hamilto- 
nian TC, neglecting the off-diagonal interaction terms and eliminating the diagonal one-body 
terms, one arrives at the Bose-Hubbard Hamiltonian (OQ). The one-body Hamiltonian reads 
h(x) = \ + V{ x ), where V(x) is a symmetric double- well potential. The left- and right- 
localized Wannier functions are obtained as linear combinations of the ground (gerade) and 
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first-exited (ungerade) states of h(x), i.e., 



<t>g(x) + <j> u (x) 

V2 



and 



= and 

oo 
oo 



are real- valued functions. The Bose-Hubbard parameters are given by U = A /_ dx(j)\{x) 
and J = — dxcj)L(x)h(x)4>R{x) . 



There is an interesting symmetry connecting the Bose-Hubbard Hami 



repulsive H{U) and attractive H(—U) interactions of equal magnitude 19 
defining the unitary operator (transformation) 



tonian ([T]) with 



2(j. Explicitly, 



it is simple to see that [19|, |20|] 



R = {b L ^ b L ,b R ^ -b R j , 



RH(U)R= -H(-U). 



(2) 



(3) 



on the evolution in 



What is the impact of the symmetry (j2J) and the resulting relation 
time of attractive and repulsive bosonic Josephson junctions? 

We consider a system of N bosons initially prepared as mentioned above in, say, the 
left well, [AT, 0) = "^ft(^l) \vac). Its evolution in time is simply given by e~ lH( - u ^ \N,0). 
Then, the "survival probability" of finding the bosons in the left well as a function of 



time is defined as pL(t;U) = J^dx ( N, 



HH(U)t 



-iH(U)t 



N, > , where 



the expression (...) is the system's density. Plugging into the "survival probability" the 
expansion of the field operator and after some algebra the final result reads: 



p L (t;U) 



X 



1 



N,0 



dx(p 2 L (x] 



1 



dx(p 2 L (x] 



1 

N 



cos[H(U)t]b[b L cos[H(U)t] + sm[H(U)t)b{b L sm[H(U)t] 



N,0 



(4) 



In obtaining the r.h.s. of (jlj) we made use of the facts that the expectation value of b^bi 
(hermitian operator) is real, and that f° dx<j)L(x)(f)ji(x) = 0. 

Employing the R symmetry (j2J) to the Pl(1;;U) matrix element (jlj) and making use of 
relation Q, one immediately finds that 



p L {t;-U)=p L {t;U), 



(5) 



which concludes our first proof. In other words, the "survival probability" of bosons is 
identical for attractive and repulsive interactions (of equal magnitude) within the Bose- 
Hubbard model. We emphasize that the result (jSJ) holds at all times t. 
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Next, we discuss the impact of the symmetry (|2J) on the eigenvalues of the reduced 
one-body density matrix within the two-site Bose-Hubbard model. The eigenvalues of the 
reduced one-body density matrix of a Bose system determine the extent to which the system 



is condensed or fragmented 



21 



22] . For the two-site Bose-Hubbard problem the reduced 



one-body density matrix can be written as a two-by-two matrix: 

\Pl R (t;U) P R R (t;U)J 

where p LL {t; U) = (V, e +iA( - u H j L b L e~ ifl( - u ^ N,0j, and p LR {t; U) and p RR {t; U) are given 
analogously. Plugging the symmetry (T5]) into each of the matrix elements of p^(t; U) and 
making use of relation one can straightforwardly express the reduced one-body density 
matrix for attractive interaction as follows: 

\-PLR(t;U) p RR (t;U) j 

Obviously, the matrices ([6]) and ([7j) have the same characteristic equation, and hence the 
same eigenvalues. We have thus shown that, within the Bose-Hubbard model, the eigenvalues 
of the reduced one-body density matrix do not depend on the sign of interparticle interaction, 
which constitutes our second proof. Again, this result holds for any time t. In particular, the 
Bose-Hubbard model attributes identical condensation and fragmentation levels to attractive 
and repulsive systems. 

To illustrate the above findings we plot in Figs. [1] and [2] the "survival probability" and 
occupation numbers, respectively, as a function of time for a repulsive and attractive Bose- 
Hubbard dynamics with the ratio of parameters ^j- = 0.226. These parameters draw from a 
double-well potential V(x) formed by connecting two harmonic potentials V±(x) = |(x±2) 2 
with a cubic spline in the region \x\ < 0.5, and for |Ao| = 0.0129. The number of bosons 
is iV = 20. The Rabi oscillation period is t Ra u = n/J = 140.66. As expected, based on 
Eq. (jSJ) and Eqs. fl6|7j) . the Bose-Hubbard dynamics for attractive and repulsive junctions 
are identical. 

We now move on to the dynamics computed with the full many-body Hamiltonian Ji. 
Recently, we have reported on the numerically-exact solution of a one-dimensional repulsive 
bosonic Josephson junction, which has allowed us to unveil novel phenomena associated 
with the quick loss of the junction's coherence fisl . We use here the same method to 



compute the time evolution of the system with the full many-body Hamiltonian 7i. We solve 
the time-dependent many-boson Schrodinger equation by using the multiconfigurational 



time-dependent Hartree for bosons (MCTDHB) method [23( (also see [2J, |25j), in which 
a novel mapping of the many-boson configuration space in combination with a parallel 
implementation of MCTDHB are exploited 26j|. This allows us to report, among others, on 
the first numerically-exact results in literature of a bosonic Josephson junction for attractive 
interaction, thus matching our recent calculations on repulsive bosonic Josephson junctions 
fis| . The results of the computations with the full many-body Hamiltonian are collected in 
Figs. [Hand EJ It is clearly seen that the dynamics of the attractive and repulsive junctions 
are distinct from each other, and that each is different from the Bose-Hubbard dynamics. 

Let us analyze these findings. We first note, in the context of the above analytical re- 
sults on the Bose-Hubbard dynamics, that the full Hamiltonian 7i does not possesses the 
symmetry (j2J) connecting the dynamics of attractive and repulsive systems. This is because 
Ji contains off-diagonal interaction terms as well as all other terms neglected in the Bose- 
Hubbard Hamiltonian ([IT) . From this "mathematical" discussion alone, we do not expect the 
dynamics of attractive and repulsive junctions to be equivalent as found above for the Bose- 
Hubbard dynamics. What do we expect on physical grounds? Intuitively, we know that 
attractive bosons like to be together, whereas repulsive bosons lean to separate from one an- 
other. These tendencies are exactly what we see in Fig. HJ The full-Hamiltonian's "survival 
probability" is larger (smaller) for attractive (repulsive) interaction than the Bose-Hubbard 
"survival probability", at least up to tjt^ a u = 1.5. In other words, the Bose-Hubbard "sur- 
vival probability" underestimates the "survival probability" for attractive and overestimates 
it for repulsive interaction for short and intermediate times. For longer times, as seen in 
Fig. (U the dynamics becomes more complex and anticipating the differences between the 
exact and the Bose-Hubbard dynamics cannot rest on the above-mentioned physical intu- 
ition alone. Finally, Fig. [2] presents a complementary picture of the dynamics of occupation 
numbers. It has been shown in 18) that the Bose-Hubbard dynamics underestimates frag- 



mentation and overestimates coherence of repulsive bosonic Josephson junctions. We may 
analogously anticipate that the reverse happens with attractive interactions, which indeed 
is the physical picture unveiled in Fig. [2j 

Let us conclude. We have shown, both analytically and numerically, that a symmetry 
present in the two-site Bose-Hubbard Hamiltonian dictates an equivalence between the evo- 
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lution in time of attractive and repulsive bosonic Josephson junctions. The full many-body 
Hamiltonian does not possess this symmetry and consequently the dynamics of the attractive 
and repulsive junctions are distinct. The Bose-Hubbard dynamics underestimates the "sur- 
vival probability" and overestimates fragmentation of attractive one-dimensional bosonic 
Josephson junctions and the reverse is true for repulsive junctions. Note that the param- 
eters used here are within the range of expected validity of the Bose-Hubbard model for 
Josephson junctions [4]. The clear deviations from the numerically-exact results show that 
criteria for the validity of the Bose-Hubbard model which have been derived for static junc- 
tions cannot be transferred for dynamically evolving junctions (also see jl8]). The present 
investigation of attractive versus repulsive junctions sheds additional light on the restrictions 
of the Bose-Hubbard model to describe dynamics. 

As an outlook, we mention that an analogous symmetry to (j2J) can be found for the 
multi-site Bose-Hubbard model. Consider the multi-site one-dimensional Bose-Hubbard 
model: H BH (U) = -J^^i+i + f E, tyfykh- Then > RbhH bh (U)R bh = -H BH (-U) 
where Rbh = |&2j — > &2j,&2j-i — ► — fyy-ij- The extension to the Bose-Hubbard model of 
orthorhombic lattices in higher dimensions is straightforward. It would be interesting to 
search for the consequences of this symmetry in the dynamics of attractive and repulsive 
bosons in a lattice. 
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FIG. 1: (Color online) Bose-Hubbard versus full-Hamiltonian, numerically-exact dynamics of at- 
tractive and repulsive Josephson junctions. Shown is the "survival probability" as a function of 
time, pi,(t), computed with the full many-body Hamiltonian for attractive [solid gray (green) line] 
and repulsive [solid black (blue) line] interaction. The Bose-Hubbard result [dashed (magenta) line] 
is for both attractive and repulsive interactions. The parameters used are N = 20, ^j- = 0.226, 
|Ao| = 0.0129, and t/j a fej = 140.66. The inset shows the ratio p" L /p^ of attractive to repulsive 
"survival probabilities" as a function of time. The black-solid line is the full-Hamiltonian results 
which exhibit a complex dynamics, whereas the dashed-magenta line is the Bose-Hubbard result, 
showing no dynamics at all. All quantities are dimensionless. 




8 




12 3 

tn Rabi 



FIG. 2: (Color online) Bose-Hubbard versus full-Hamiltonian, numerically-exact dynamics of at- 
tractive and repulsive Josephson junctions. Shown are the occupation numbers of the reduced 
one-body density matrix as a function of time, n\ (t), computed with the full many-body Hamil- 
tonian for attractive [solid gray (green) line] and repulsive [solid black (blue) line] interaction. The 
Bose-Hubbard result [dashed (magenta) line] is for both attractive and repulsive interactions. The 
two-site Bose-Hubbard dynamics has two occupation numbers only. The full-Hamiltonian dynam- 
ics has many occupation numbers. It is seen that the occupation numbers n^ 2 (t) are essentially 
zero. The parameters used are the same as in Fig. [TJ All quantities are dimensionless. 
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